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Abstract: Spin-2 particles and in particular gravitons are predicted in many new physics sce- 
narios at the TeV scale. Depending on the details of models such new states might show up 
as a continuum, massless particles, or TeV scale resonances. Correspondingly, very different 
discovery signatures should be exploited, from the search of excesses in events with multi jets 
and large missing transverse energy, to resonances in weak boson or jet pair productions. We 
present a very general and flexible implementation in MadGraph/MadEvent of spin-2 par- 
ticles interacting with the standard model particles via the energy momentum tensor, which 
encompasses all of the most popular TeV scale models featuring gravitons. By merging matrix 
elements with parton shower, we can generate inclusive samples of graviton+jets at the hadron 
colliders in several scenarios (ADD, zero-mass graviton and RS). We compare and validate our 
results against the corresponding next-to-leading order QCD calculations. 
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1 Introduction 

Currently, the Standard Model (SM) of particle physics matches to a great accuracy a plethora of 
the experimental data. However, there are many reasons to expect new physics beyond the SM 
appearing at the LHC (TeV) energy scale, such as the demand for a dark matter candidate or the 
problem on the large hierarchy between the electroweak and Planck scales. Many models attempt 
to provide a solution to these and other issues by extending the SM to larger symmetries (such 
as in Supersymmetric or Little Higgs models) or larger space-time (extra space dimensional 
models), with the condition that they reproduce the SM at low energies, while new physics 
effects might only become apparent at high energies. The mechanism of how these new degrees 
of freedom appear at the TeV scales is particularly intriguing for extra dimensional theories. 
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Take the flat extra dimensional models (ADD) [1-3] for example. In this D = 4 + 6 dimen- 
sional model, only gravitons (G) can propagate in the D dimensional space-time, while all the 
SM particles are confined onto the usual 3-|-l-dimensional space-time. Assuming for simplicity 
that the additional (5-dimensional space to be compactified on the torus with a common radius 
R, the Planck scale (Mpi) can then be related to the fundamental scale Ms as follows: 

M|i = 8ttR^M^+'^ . (1.1) 

If one chooses a large compactification radius {R ^ 1/Mpi), the hierarchy between the Planck 
and the electroweak scale can be understood by the fact that the fundamental scale is now 
allowed to be closer to the weak scale, i.e., Mg ~ TeV, where the effects of extra dimensions 
will start manifest themselves through an apparent non-conservation of energy and momentum 
in the usual 3-1-1 dimensions. 

In addition to the ADD model, other TeV scale theories on gravity have been proposed 
that attempt to explain the hierarchy problem in similar or different ways. For instance, if 
a huge number of particles are in a hidden sector that interacts only gravitationally with the 
SM particles, the Newton's constant might be affected and the fundamental gravity's scale 
could then be in the TeV region (TeV scale Massless Graviton Model or MGM for short) [4-9]. 
Alternatively, in warped extra dimensional models proposed by Randall and Sundrum (RS) in 
1999 [10, 11], the large hierarchy is addressed by the introduction of a geometrical exponential 
factor in the metric, a solution that does not need the extra dimension to be significantly larger 
than the Planck length. 

The search for this kind of new physics at the TeV scale is today among the major tasks 
of the LHC. The most natural signatures involve missing energy in association with jet(s) {e.g., 
in the ADD and MGM models), or spin-2 resonances in lepton or jet pair productions {e.g., 
in the RS model). The experimental analyses of such signatures demand accurate simulations 
within Monte Carlo (MC) event generators not only to design efficient signal selection strategies 
but also to extract information on the new physics parameters (masses, coupling strengths, 
compositeness, structures and spin). Several tools exist nowadays that are able to perform such 
simulations for a wide class of models [12]. In this work, we will focus on MadGraph [IS- 
IS], where a spin-2 particle implementation at the matrix element level has been presented 
recently [16]. In this work we have extended it to full-fledged automatic event generator and 
embedded it in the MadGraph/MadEvent (MG/ME) package. This now includes a wide 
class of models involving spin-2 particles, such as: 

• weakly coupled (Planck scale suppressed) massive graviton models, with infinite number of 
nearly continuously distributed Kaluza-Klein (KK) modes of the graviton and the overall 
integrated effects are significant at the TeV scale, e.g., the ADD model; 

• TeV scale massless graviton models, in which the graviton effects are negligible below 
the fundamental scale of quantum gravity, i.e., the TeV scale, due to the running of the 
gravitational constant, e.g., in theories with many particles that interacts with the SM 
particles only gravitationally; 
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• strongly coupled (TeV scale) massive graviton models, in which the graviton KK modes 
are widely separated and the lowest one gets mass at the TeV scale and can decay into 
the SM particles within detectors, e.g., the RS model. 

To show and validate the capabilities of the new implementation we have chosen to study the 
most important signatures in representative models of the three classes above. For the ADD and 
MGM models, the gravitons couple weakly to the SM particles and thus can appear as missing 
energy signals at colliders. The main searching channel is pp{p) — s- G + jet(s), in which at least 
one accompanying hard jet is needed. For the RS model, resonance production pp{p) — )■ G is 
expected, with G decaying into e^e~ or /U+/i~ for example. In fact, for completeness we further 
consider the G plus mono-jet channel to compare the RS model with the ADD and MGM 
results. The corresponding NLO QCD calculations to those channels have been performed in 
Refs. [17-19]. 

In this paper we have extended the NLO QCD calculations in [17] to the cases of the 
RS and MGM models. While current NLO fixed-order results are useful to reliably predict 
cross sections and observables involving at most one jet, multi-jet based observables and more 
exclusive quantities in general can have an important role in the analyses. For example, NLO 
QCD calculations for G + mono-jet production can only provide LO distributions for the 2nd 
jet G + jet. Merging matrix elements with a parton shower can, on the other hand, reliably 
predict multi-jet final states and generate samples that can be directly used in experimental 
simulations. Several methods of merging or matching schemes are now available [20-22] which 
have been shown to be in good agreement [23]. 

Here, we present the first results from the parton shower merging with matrix elements, 
for spin-2 graviton productions at hadron colliders, where the new implementation of graviton 
in MadGraph/MadEvent and A;_l— MLM merging scheme [24] are used. We further present 
a detailed comparison between the matching and the relevant NLO/LO results [17-19]. We 
also include the results for an irreducible and representative background process from hadro- 
productions of Z obtained with MCFM [25]. 

The article is organized as follows. In section 2 we review briefly the three models men- 
tioned above. The complete implementation of spin-2 particle in MadGraph/MadEvent is 
detailed in Section 3. In Section 4, we give an introduction for parton shower matching in Mad- 
Graph/MadEvent. The numerical results of the matching and the NLO/LO calculations for 
G hadro- productions are presented in Section 5. We draw our conclusions in Section 6. 

2 Description of the models 

2.1 Extra space Dimensions: ADD and RS models 

Various extra dimension models which have been proposed so far can be divided into two major 
classes according to the geometry of the background space-time manifold. The first one includes 
the ADD and its variants, which extend the dimension of the total space-time to D = 4-1-5, with 
a factorisable metric and large size of the compact extra dimensions l/Mpj). The second 
one includes the 5-dimensional RS model [10, 11] and its variants, in which a warped metric is 
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introduced along the 5-th dimension the size of the extra dimension can be at the order of the 
Planck length. 

In both classes of extra dimension models Kaluza-Klein (KK) towers of massive spin-2 gravi- 
tons appear, that can interact with the Standard Model (SM) fields. The effective interaction 
Lagrangian is given by [26, 27]: 

Ant = -^E^^"^'"V' (2.1) 

n 

where T^'^^^'^ is the n-th graviton KK modes, and A is the relevant coupling scale. In the ADD 
model 

A = Mpi = Mpi/\/8^ ~ 2.4 X 10^*^ GeV, (2.2) 

where Mpi is the 4-dimensional reduced Planck mass, while in the RS model 

A = e-^^=Wpi , (2.3) 

where A: is a scale of order of the Planck scale and rc is the compactification radius. By tuning 
krc one can place A at the electroweak scale. 

In Eq. (2.1), is the energy-momentum tensor of the SM fields, 



(2.4) 



where g^'^ is the metric and j/^*^ = diag(l, —1, —1, —1) is the Minkowski value. Note the difference 
on the effective strength of the gravity/SM coupling for each class of models. 

In the ADD model the individual KK resonances have masses equal to m(^) = \n\/R. The 
mass gap between neighboring modes Am = is thus small for 6 not too large. Quantitatively 
one finds Am 20 keV, 7 MeV and 0.1 GeV for Ms = 1 TeV and <5 =4, 6 and 8, respectively [27]. 
The discrete mass spectrum can be approximated by a continuum with a density of states 
diV = p{m)dm [27, 28], where 



Mpi , „ 27r^/2 



p{m) = Ss^ij^^m^ \ andSs^i = ^jjj^. (2.5) 



In other words, the 4-dimensional graviton appears as an infinite sum of excited states of the 
graviton. Current terrestrial test of gravity set a limit on Mg > 3.6 TeV for 5 = 2 [29]. Further 
constraints have been derived from astrophysics and cosmology, in particular for 6 < 4. However, 
they can be evaded in specific models [30-32] and do not lessen the importance of collider 
searches for extra dimensions. At high energy colliders, both virtual graviton exchange between 
SM particles and real graviton emission provide viable signatures of large extra dimensions. 
Since the coupling of gravitons with matter is suppressed oc 1/Mpi, direct graviton production 
gives rise to missing energy signals. Searches for the ADD graviton production have been 
performed in the processes e+e~ — )• 'yiZ) + E"^^^^ at LEP and pp — )• 7(jet) -|- p^^^^ at the 
Tevatron. The combined LEP limits [33] read Ms > 1.60, 1.20, 0.94, 0.77, 0.66 TeV, for 6 = 
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2,- • • ,6 respectively, while Tevatron searches exclude Mg > 1.40, 1.15, 1.04, 0.98, 0.94 TeV, for 
(5 = 2,- • • ,6 respectively [34-36]. 

In the RS model, the mass of the nth graviton KK excitation mode at the electroweak scale 
is given by 

mn = kxne-^^'^ = mi — , (2.6) 

Xl 

where the the rath roots of the first order Bessel function. The graviton sector of the RS 

model is completely determined by the two parameters mi and A. Current constraints [37-39] 
for the parameters of the RS model are from the theoretical requirement, the low energy precise 
measurement and also the data from Tevatron, from which 0.01 < fc/Mpi < 0.1 and A < 10 TeV. 
The relation between k/M-p\ and A is given by: 

A=!^. (2.7) 

Xl k 

At variance with ADD, the lightest RS massive graviton can have a mass of TeV scale, and may 
be produced copiously at the LHC. More importantly, it has much larger couplings to the SM 
particles than the ones in the ADD model, which allows the graviton to decay into observable 
particles and hence be detected as a resonance. 

2.2 The Massless Graviton Model 

In the previous section, we have seen that if Nature entails more than 4 space-time dimensions, 
gravity could become strong at much lower scale than the Planck scale. However, it has recently 
been shown that even for a 4-dimensional model, TeV scale gravity effects could appear [5, 8]. 

For instance, Ref. [5] suggests that strong gravity should emerge at the TeV scale if there 
are huge number of hidden sector particles that interact only gravitationally with the Standard 
Model particles. In Ref. [8], it has been suggested that the emergence of TeV scale quantum 
gravity may be interpreted as renormalization of the effective gravity coupling due to hidden 
sector particles, based on the one- loop computation [40]. Common to these models [4-9] is the 
existence of order {M-pi/TeV)'^ ~ 10'^^ hidden particles at or below the TeV scale, that can only 
be probed by the gravity interactions. 

The effective Lagrangian for this type of theory can be expressed by the effective Fierz-Pauli 
Lagrangian combined with the interaction term of the graviton to the Standard Model particles: 

where T^^'^ is the massless graviton and T^'^ is the SM energy-momentum tensor given by 
Eq. (2.4). The graviton coupling has the same form as that of Einstein's theory of general 
relativity, except for the effective Planck mass M(/x^) which is assumed to be in the TeV range 
at energy scale /Xj, ~ TeV. One can note the similarity between Eq. (2.8) and the Eq. (2.1). As a 
consequence, similar Feynman rules are expected. A detailed set of Feynman rules for the extra 
dimensional theories can be found in Refs. [26, 27], and for the MGM Model in Ref. [41]. 

Presently, the MGM model is less constrained than ADD or RS models. It has been shown 
that the cosmic rays experiment AGASA sets a bound of 550 GeV on the effective Planck mass 
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M{fj,^) in 4-dimensions [42]. However, more recently results on monojet in addition to graviton 
emission [41, 43] predict this limit to be at most 0(5 TeV) at the LHC (14 TeV). 

3 Graviton implementation into MadGraph 

Fortran subroutines to calculate helicity amplitudes with massive gravitons were added to the 
HELAS (HELicity Amplitude Subroutines) library [44, 45] recently [16]. Thanks to these new 
subroutines which encode the Lorentz structure of the new couplings, the ADD and RS mod- 
els could be implemented in MadGraph after a few modifications. The implementation in 
MadGraph/MadEvent, however, was not complete due to a few technical limitations. 

First, MadGraph version 4 can only generate Feynman diagrams with up to 4-point ver- 
tices, the diagrams involving the 5-point vertex (4 gluons-graviton) have therefore to be added by 
hand. Although it was fine for previous studies on graviton production with up to 2 jets [17, 46], 
it is inconvenient for graviton production with higher multiplicity jets. 

Secondly, gravitons are densely distributed in the ADD model (see Eq. (2.5)) and a special 
treatment for the propagator is called for. 

Third, for the MGM model it is necessary to have the off-shell graviton current subroutines 
modified due to differences on the graviton propagator and polarization summation for the 
massless and massive cases. 

We present hereby the details on how each of the above issues has been tackled to achieve 
spin-2 graviton generation within MadGraph/MadEvent in a full automatic way. 

3.1 Five-point vertices 

The interactions (2.1) and (2.8) entail a five-point vertex, ggggT: 



with the coupling constants GT = — and GC = g^- In HELAS such a vertex can be encoded as 
GGGGTX computing the portion of the amplitude of the ggggT vertex from four Gluon polarization 
vectors and a Tensor boson wavefunction corresponding to the color structure J«*<=ycde^ 
figure 1 the diagram on the left-hand side. In order to implement this vertex in MadGraph we 
reinterpret it as scattering in the t, u and s channels of non-propagating octet tensor boson 
tA as an auxiliary particle. With we can reduce the portion of the 5-point vertex with 
color structure into 3-point vertices: ggtA and tAtAT, which can then be generated by 

MadGraph automatically, figure 1 the diagram on the righ-hand side. Note we assign a fiow 
to tA, i-G., we treat tA different from its antiparticle in MadGraph , and forbid the interaction 
with the antiparticle to avoid appearance of additional diagrams for gg — t- gg. 

After introducing particle tA, we add several new HELAS subroutines as specified on table 
1, by matching the amplitude of the right side in figure 1 with the portion of the left side of the 
same color structure, i.e., Eq. (A. 3). The interested reader will find a detailed description in 
Appendix A. 

The new HELAS subroutines have been tested by using the QCD gauge invariance and 
the general coordinate transformation invariance of the helicity amplitudes. Finally, we mention 
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Figure 1. Reduction of tlic portion of tlic 5 point vertex with color term /"■'"^/'^"■'^into 3-point vertices 
by introducing auxiliary tensor particle Ia- T represents the graviton. 



Sub-routine 


Summary 


VVTCXX 


computes the amplitude of the ggt^ vertex. 


JVTCXX 


computes the off-shell vector current J made from the interactions of a 
Vector gluon and an auxiliary Tensor boson by the ggi^ vertex. 


UVVCXX 


computes the off-shell tensor current U for the auxiliary tensor tA, made 
from two flowing-out Vector gluons by the ggtyi vertex. 


TTTXXX 


computes the amplitude of the tA^ vertex. 


UTTAXX 


computes the off-shell non-propagating tensor current U, made from the 
flowing-out graviton Tensor and the auxiliary Tensor by the Ia ^aT ver- 
tex. 


UTTBXX 


computes the off-shell graviton tensor current U, made from the two 
flowing-out auxiliary Tensors by the Ia t^T vertex. 



Table 1. Summary of new vertices introduced for the implementation of the 5-point vertex in the 
HELAS structure. 



that the above technique has also been applied for a new implementation of the Higgs effective 
theory (HEFT) in MadGraph, see details in Appendix D. 

3.2 Summation over the graviton mass spectrum 

The 4-dimensional graviton can be approximated as an infinite tower of KK graviton modes 
within the ADD theory. A MC simulation has therefore to take into account an extra phase- 
space integration to perform: the amplitude integration over the graviton mass density function 
(Eq. (2.5)). 

For virtual graviton exchange, the tower of KK graviton modes leads to the summation of 
their propagators, which is ultraviolet (UV) divergent since ADD is only an effective theory. 
One usually introduces an UV cutoff for the highest KK modes and replaces the summation by 
factors in some effective scale [26]. This can be easily implemented in MadGraph by setting 
the graviton mass to zero and the extra dimension scale to be the effective one. 

Nevertheless, it is more complicated for the real graviton production case, for which we 
need to actually perform the integration described above. The integration is carried out by 
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Sub-routine 


Summary 


PXXXXX 


stores the helicity and momentum of the auxihary particles xi and X2- 


TPSXXX 


computes the ampHtude of xiX2'^ vertex from a set of parameters spec- 
ified by the theory. 


UPSXXX 


computes an off-sheh graviton tensor current U, made from the flowing- 
out auxiliary Scalar and Pesudo-particles, by the xi x-^l vertex. 



Table 2. Summary of new vertices introduced for the implementation of the mass integration in Mad- 
Graph/MadEvent within the framework of the ADD model. 



the introduction of two massless colorless auxiliary particles: a scalar xi, and a pseudo-tensor 
X2- MadGraph is then modified to include a new particle type P, which is only involved in 
the interaction with graviton: x\ X2Y.. Instead of generating process such as — t- T + X in 
MadGraph, one can generate — t- x\X2 + X. The phase spaces are related by 

d$(xiX2+X) = — ^d«>(r + X). (3.2) 

Using Eq. (3.2), we can perform the mass integration with the help of the new HELAS subrou- 
tines for the x\ X2I vertex summarized in table 2. Details on each new sub-routine are reported 
in the Appendix B. 

To conclude we mention that the new HELAS subroutines and the new implementation 
of the propagators have been tested by comparing with previous implementation [16] always 
finding excellent agreement. 

3.3 Implementation of the massless graviton model 

A well-known property of spin-2 particles in the framework of Einstein's general relativity is the 
so called van Dam-Veltman discontinuity [47], which states that the massless propagator for a 
graviton cannot be obtained by the massive one through a smooth limit. Such discontinuity is 
due to the fact that the massive spin-2 particles have five states of polarization, while massless 
ones have only two of them. This can be observed on part of the the massive/massless graviton 
propagators [26, 4f]: 

A=±2,±l,0 

12 2 

'n^l|3PvPa + rjuaPt^Pp) + ^ {v^^ly -I 2Pi^P^)^^l(^P jPo'Pp)^ (3-3) 



B'^^^'r^ E e/..(P,A)e„;3(p,A) 

A=±2 
1 



), (3.4) 
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where = {p^,p) is the four-momentum of the graviton, fj^,^ = r]^y — ^^^^ ^^^^ and = 
(jfi , —p). One can always check that the graviton propagators satisfy B^y^^iiP^ = 0, and lead to 
the correct number of spin degree of freedom by computing the traces: 

Observe from Eqs. (3.3) and (3.4) that B^^^J^'^ in the limit when rriT — t- is different than 

omassless 

This property has to be taken into account for the HELAS routines that includes massless 
gravitons. The corresponding modifications have been introduced for: i) the tensor wave function 
subroutine and ii) all the subroutines with an off-shell graviton. 

The new subroutines have been tested through gauge invariance, and by comparing results 
of total cross sections for various 2—7-2 processes against the massive graviton case. The 
reason why we are allowed to perform the comparison is that although the discontinuity on the 
graviton propagator exists, it does not occur at the cross section level for processes relevant at 
colliders [28, 41]. 

Moreover, the implementation of the running of the Planck mass has been done in Mad- 
Graph/MadEvent during the generation of the events: we have applied a minimum cut of 
H™™ > M{p,^) on the sum of the jets transverse momentum. In other words, a Heaviside step 
function-like behavior is added to the cross section, i.e., for collisions with s, (— t), {—u) < 
M{ni,), gravity contributions are so weak that the cross section can be approximated to zero. 

4 Parton shower merging with tree level matrix elements 

As the parton center of mass energy increases, hadron collision events with large jet multiplicities 
become more probable. As a result, accurate simulations require to correctly account for the 
presence of QCD radiation, which might modify or alter the leading order predictions for the 
relevant observables. 

In parton shower MC programs like Pythia [48] and Herwig [49] additional jets are usually 
obtained in the collinear and soft approximation. Hard and widely separated jets are therefore 
poorly described in this approach. On the other hand, tree-level fixed order calculations can 
provide reliable predictions in the hard region, while failing in the collinear and soft limits. 

To combine both descriptions and avoid double counting or gaps between samples with 
different multiplicity, an appropriate matching method is required. Several algorithms have 
been proposed over the years: the CKKW method, based on shower veto and therefore on event 
re- weighting [20, 21] and the MLM-based scheme, based on event rejection [23]. 

For SM processes, in particular for the weak vector boson hadro-productions, matching 
schemes have been extensively applied and compared to the available data from the Tevatron [50]. 
Furthermore, studies have been performed for new physics in order express the significance of 
the matching, such as the notable reduction of the dependence on parton shower parameters [24]. 

In this study, we use the A;_|_-MLM matching scheme implemented in Mad Graph /MadE vent 
and interfaced to Pythia for parton shower and hadronization. Detailed comparisons of sev- 
eral event generators that include matching schemes have been performed in Refs.[12, 23, 51], 
showing a good overall agreement. 
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In the k± matching scheme, matrix element multi-parton events are produced with a mini- 
mum separation k± cutoff of Q^^- For every event, the final-state partons are clustered accord- 
ing to the k± algorithm. Only clusterings corresponding to the Feynman diagrams provided by 
MadGraph are allowed. In order to closely mimic the behavior of thin the RS model, massive 
grae parton shower, the k± value for each clustering vertex corresponding to a QCD emission 
is used as renormalization scale for in that vertex. For the central hard 2 — >• 1 or 2 — )• 2 
process, the transverse mass = p'± + fn^ of the particle(s) produced in the central process is 
used as a factorization scale, as well as a renormalization scale. 

Subsequently, this event is passed to the parton shower MC simulator. Before hadronizing 
or decaying, the final partons are clustered into jets using the k± algorithm with a jet cutoff of 
Qmin Qmin- '^^^ J^^^ then Compared to partons. They are considered to be matched if 
A;_|_ (parton, jet) is smaller than the cutoff Q^^i^- An impossibility of matching all partons with 
jets, results on the rejection of the event. For events with parton multiplicity smaller than the 
highest multiplicity, the number of jets must be equal to the number of partons. 

4.1 Choice of matching parameters and cuts 

As mentioned above, the aim of this paper is to study inclusive graviton production, i.e., pp (p) — )■ 
G + n jet(s) at the LHC and Tevatron, within the three models: ADD, MGM and RS. 

For the ADD and MGM models, missing energy with hard jets will be the signal at hadron 
colliders to identify real emissions of gravitons. The main background is pp (p) ^ Z + {n-)j 
with Z boson decays into neutrino pairs. In the following, we consider matrix elements with 
parton multiplicity n from 1 to 3 for the matching with parton shower for the signal processes. 

To suppress SM backgrounds in the LHC graviton searches [52], we require 

P|P'"" > 500 GeV . (4.1) 

Jets are defined by the kj- algorithm, with the resolution parameter set to D = 0.6, and are 
required to satisfy |?7j| < 4.5 and > 50 GeV. 

Regarding the matching parameter at the LHC, the separation cutoff is set as: 

Q^l>A5GeY, (4.2) 
gj;:L>50GeV. (4.3) 

At the Tevatron analysis we use the same settings as in the recent CDF study [34, 36], i.e., 

P^'^^ > 120 GeV , > 150 GeV , and \rij\ < 1 . (4.4) 

Here, jets are defined by the kj- algorithm with D = 0.7, and are required to satisfy \r]j\ < 3.6 
and Pj, > 20 GeV. A second jet with Pt > 60 GeV is vetoed. The separation cutoffs used for 
the matched analyses are: 

Q^l>20GeY, (4.5) 
gi^i„>30GeV. (4.6) 

In the RS model, massive gravitons can decay quickly into visible particles within the 
detector. Thus we first generate a full inclusive sample for G productions at the hadron colliders. 
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for which we take into account the matrix elements with parton multiphcity n from to 2. We 
then let G in the inclusive sample decay into e'^e~ or ^~ , keeping spin correlations exact. 

To allow a thorough comparison with the NLO QCD calculations on graviton mono-jet 
productions of Ref. [17], and in analogy with the analysis of the ADD and MGM models, we 
generate a semi-inclusive sample for G-jet productions in the RS model with the same cuts and 
settings as for the ADD and MGM models. For simplicity, we do not ask for G decay here, as 
the aim is to validate the matching method against NLO QCD calculations. 

5 Results 

Let us fix now the model parameters. Considering the up-to-date experimental constrains as 
mentioned in Section 2.1, for the ADD model, we choose A = 5 TeV for the LHC, A = 1 TeV 
for the Tevatron, with = 2, 4, 6. For the RS model, we take nii = 1 TeV and 100 GeV, with 
A = 3 TeV. For the MGM model, we set M(/x^) ~ /ij, = 1 TeV and 2 TeV. 

In the event generation with MadGraph/MadEvent , CTEQ6L1 PDF [53] is employed 
which also fixes the values of as at the Z mass, while the choices on renormalization and 
factorization scales are defined to be the transverse mass of the graviton as explained in Sec. 4. 
Correspondingly, for the LO/NLO calculations, CTEQ6L1/6M PDFs are employed together 
with the corresponding values for the strong coupling as- The renormalization and factorization 
scales are set to the transverse momentum of the graviton Pj?, and we will also show the 
uncertainties by varying the scales in the range between /2 and 2P^ . 

One relevant remark concerns the range of validity of the effective theories considered here. 
As they all assume a linearized Einstein gravity, they are valid only when the scales involved in 
the hard scattering process do not exceed the fundamental scale Ms, in which case a quantum 
gravity description is needed. 

In the following we present the matrix element-parton shower matched results and the 
LO/NLO ones, for the differential distributions: 

• Pseudo-rapidity of the leading jet (rf )] 

• Pseudo-rapidity (rp) of the graviton; 

• Transverse missing Pt for ADD and MGM models, or graviton Pt for RS model; 

• Leading and Second jet Pt; 

• Ht of the jets, which is defined as the sum of all the jet Pt in each event: 




(5.1) 



Model 


K-factor Normalization factor 


mi = 1 TeV 
RS mi = 500 GeV 


1.65 1.73 
1.75 1.82 



Table 3. K-factor (J {n lo) I <^ (lo) and normalization factor Oi^j^^o)! '^{mlm) f^i' inclusive G production in 
the RS model at the LHC, from pp ^ G NLO QCD calculations and parton shower matching with the 
matrix elements oi pp ^ G + n jet(s) with n = 0, 1, 2, A = 3 TeV, mi = 0.5 and 1 TeV for y/s = 7 and 
14 TeV LHC, respectively. 



Model 


Normalization factor 
LHC Tevatron 


6 = 2 
ADD (5 = 4 
<5 = 6 


2.05 2.40 
2.34 2.47 
2.49 1.92 


Jlifi^) = 1 TeV 
MGM Wifi^) = 2 TeV 


1.56 
1.64 


mi = 1 TeV 
RS mi = 100 GeV 


1.99 1.95 
1.81 1.73 



Table 4. Normalization factor cr(NLO)l '^(mlm) foi' semi-inclusive G production in different models at 
the LHC and Tevatron. 

5.1 Inclusive sample for the RS graviton productions 

Let us start the discussion with the results for the RS model where we first consider a full inclusive 
graviton production with graviton decay. The NLO QCD corrections for dilepton production 
in RS model via graviton production have been calculated in [18, 19] and compared to the LO 
expectation. Table 3 collects the results in terms of rather important K-factors (as given by 
^ {N LO) / (LO)) which have therefore to be included to achieve a satisfactory normalization of 
our fully inclusive samples, figures 2 and 3. 

In figure 2, fc_L-MLM matched results for pp ^ G + n jet(s) with n = 1,2 are presented. 
The graviton is decayed subsequently into a pair of leptons {i.e., either e^e~ or fi^fi") and the 
corresponding results are shown in figure 3. Here normalization factor between the A:_|_-MLM 
matched and the NLO expected cross section are also computed. These results are shown on the 
third column of table 3. Observe that the values for the NLO/LO normalization factor are very 
close to the ones found for the NLO/MLM normalization factor. In fact, the cross section after 
matching analysis is expected to be similar to the one computed for pp ^ G hy LO calculation. 
That is because LO calculation already includes the production of extra radiation through the 
PDFs (parton distribution functions). However, as the matching calculation considers explicitly 
the production of extra jets in the final state, it modifies the behavior of the distributions 
being analyzed such as the missing Pt, pseudo-rapidity and Ht distribution. As a result, we 
expect that matched sample, normalized to NLO calculation, will yield better predictions for 
observables which involve extra radiation then a pure parton shower approach. 
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Figure 2. fcj^-MLM matching results for a full inclusive sample of the RS model with A = 3 TeV. The 
matched results are given by the red and blue curves for mi = 1 TeV at the LHC with y/s ~ 14 TeV, 
and mi = 500 GeV at the LHC with -y/s = 7 TcV, respectively. No cuts were applied in these samples. 
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Figure 3. fcj^-MLM matching results for a full inclusive sample of the RS model with A = 3 TeV and 
the graviton decayed into a pair of leptons. The red and blue distributions are related to the mass of the 
first mode of the graviton and the collision energy. The red curves are for rrii ~ 1 TeV at the LHC with 
y/s ~ 14 TeV, and the blue ones for mi = 500 GeV at the LHC with -^s — 7 TeV. No cuts were applied 
in these samples. 

5.2 Semi-inclusive sample for graviton-monojet studies 

Figures 4-9 show all the differential distributions for the comparison between the NLO and 
k±-MLM matched results. The NLO uncertainty bands are also presented by setting the renor- 
malization and factorization scales to fip = ytR = /io/2, /^o and 2^o with /iq = ■ The A;_l-MLM 
matched results are normalized by the NLO ones with /i_p = /Ur = fJ-o, as shows table 4. 

For ADD and MGM classes of models it is much more interesting to obtain a comparison 
between G + jet at NLO and G + n jet(s) with the A;_l-MLM matching for n = 1, 2, 3. In order 
to compare all the three classes of models, we also present the same analysis for the RS model. 
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Figure 4. NLO/fc_L-MLM matching comparison for the ADD model at the LHC {y/s = 14 TeV) with 
A = 5 TeV. The NLO results are given by the red, green and purple bands for S = 2,4,6, respectively. 
The matched results are given by the black curves with the same parameters as for the NLO ones. The 
dominant Z — > i^i/ background is also shown as a reference by the blue curve. The fcj^-MLM matched 
curves are normalized by the NLO results, and the normalization factors can be found in table 4. 
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Figure 5. NLO/fcj^-MLM matching comparison for the ADD model at the Tevatron {y/s ~ 1.96 TcV) 
with A = 1 TeV for S = 2,4, 6. The matched curve is normalized by the NLO results and the normalization 
factors can be found in table 4. 

We start by discussing the results for the ADD model. Prom figure 4, one can immediately 
see that the NLO and A;_l-MLM matched results for the LHC agree quite well on the shapes in 
pseudo-rapidity of the leading Pt jet, pseudo-rapidity of the graviton and transverse missing Pt- 
For the leading jet Pt distributions, the agreement is satisfactory for large Pt ^ 800 GeV region, 
while for Pt ~ 500 GeV, the NLO curves drop down too quickly and become unreliable. That is 
due to the cut chosen as in Eq. (4.1) which leads to inconsistency on leading jet Pt lower bound 
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Figure 6. Graviton mass distribution at the LHC (left) and Tevatron (right) for the ADD model with 
(5 = 2,4, 6, respectively. The NLO results are shown by the red, green and purple bands, and the matched 
ones are shown by the black curves. 

at LO and NLO calculations, while /c^-MLM matching gives reliable results. The same argument 
also applies for Ht distributions when Ht ~ 500 GeV. However, at large Ht ^ 1500 GeV, the 
A;_i_-MLM matching tends to give harder distributions, because the contributions from the matrix 
elements of G + 3 partons are also included. 

Figure 5 presents similar results for the Tevatron, and one can find, in general, a better 
agreement between NLO and /c^-MLM matched results. Tevatron's energy scale is much lower 
than the LHC and thus the contribution from more hard jet emission included in fc_L-MLM 
matching does not play an important role within this case. 

It is known (Sec. 3.2) that the 4-dimensional graviton is represented by a infinite sum of 
KK gravitons in the ADD model. The mass of the effective graviton is therefore given by a 
distribution described by Eq. (2.5). In figure 6 we show the results for the NLO and A:_|_-MLM 
matching comparison of the graviton mass distribution. One can easily observe the dependency 
on the number of extra dimension since the average graviton mass increases with it. Again, 
/i;_L-MLM matching results agree quite well with the NLO ones. 

The results for the MGM model at the LHC are presented in figures 7. The Tevatron 
plots are out of reach for M(/i^) = 1,2 TeV, therefore the corresponding results are not shown. 
Moreover, we apply here the cut > M(//*) as described in Sec. 3.3. One can see a kink 

or turning points in the distributions of the leading or second jet Pt and transverse missing 
Pt- The NLO and A;_|_-MLM matched results agree quite well on shapes for distributions of 
pseudo-rapidity of the leading Pt jet and graviton, transverse missing Pt and Ht of the jets. 
For the second jet Pt distribution, /c^-MLM matching gives a considerably different shape 
when compared to the NLO calculation, especially for M{^^) = 2 TeV case. For the leading 
jet Pt distribution, the NLO and A;_|_-MLM matched results for the LHC agree well for large 
Pt > M(/i^)/2, while for Pt ^ Af(/ij,)/2, the NLO curves are unreliable, again due to the 
inconsistency on leading jet Pt lower bound at LO and NLO calculations. 
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Figure 7. NLO/fc^-MLM matching comparison for the MGM model at the LHC. The NLO results are 
given by the red and green bands for M(/i^,) = 1,2 TeV, respectively. The matched results are given by 
the black curves with the same parameters as for the NLO ones. The dominant Z — > i/i? background is 
also shown as a reference by the blue curve. The fcj_-MLM matched curves are normalized by the NLO 
results, and the normalization factors can be found in table 4. 
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Figure 8. NLO/fc_L-MLM matching comparison for the RS model at the LHC with A = 3 TeV. The 
NLO results are given by the red and green bands for mi = 1 TeV and 100 GeV respectively. The 
matched results are given by the black curves with the same parameters as for the NLO ones. The 
k±-MLM matched curves are normalized by the NLO results, and the normalization factors can be found 
in table 4. 
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Figure 9. NLO/fcj^-MLM matching comparison for the RS Model at the Tevatron with A = 3 TeV, for 
mi = 100 TeV and 1 TeV, respectively. The k±-MLM matched curve is normalized by the NLO results 
and the normalization factors can be found in table 4. 

Figures 8-9 presents results for the RS model at the LHC and Tevatron similarly to the 
ADD figures showed in figures 4-5. Here, the agreement between the NLO and A;j_-MLM 
matched results is again much better at the Tevatron than the LHC for the same reasons 
pointed above. Also, fc^-MLM matching tends to give harder Ht distributions especially for 
large Ht > 1500 GeV. 
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Figure 10. NLO/fc_L-MLM matching comparison for the RS Model at the LHC with mi = 100 GeV and 
A = 3 TeV. Different cuts from the ones specified above are used to show the interesting case with large 
K-factors. 

5.3 Large K-factors 

In multi-jet final states, observables exist that can be subject to large NLO QCD corrections 
at TeV scales [54]. Examples include Hx or the Pt of the leading jet. Such observables are 
sensitive to the opening of new channels with new kinematic topologies at NLO, which are 
enhanced by parton distribution effects or by special kinematic configurations (leading to large 
logs). One is then lead to the impression that higher order corrections can be dominant and 
the perturbative series is not well-behaved. In fact, one can easily show that this is not the case 
for more inclusive quantities and this behaviour is due to the choice of the specific observable. 
In any case, the effects in some specific region of the phase space can be large and need to be 
accounted for. One method was introduced Ref. [54]. Here we show that multi-parton samples 
with /cj_-MLM matching provide a reliable description also of these kind of effects. In figure 10 
we collect some representative plots which show it. For simplicity, we consider the RS model at 
the LHC with A = 3 TeV and nii = 100 GeV, with jet cuts \r]j\ < 4.5 and > 50 GeV, and 
graviton cut > 50 GeV. Comparing the LO and NLO curves, one can see the appearance 
of very large K-factors, especially at large PjT, leading jet Pt and Ht- The matched samples 



- 21 - 



describe extremely well the NLO shapes. 
6 Summary and conclusions 

In this paper we have presented a detailed comparison between predictions from next-to-leading 
order QCD calculations and the tree level matrix-element matching with parton showers in 
the context of three representative theories featuring a graviton: ADD, RS and MGM. This 
study validates and motivates the use of a fully automatized event generator such as Mad- 
Graph/MadEvent including spin-2 particles and multi-parton jet matching for LHC simula- 
tions. Let us summarize the main aspects of our approach and the key results. 

First, we have presented the complete implementation of spin-2 particles in 
MadGraph/MadEvent event generator, which now includes: (i) new HELAS subroutines to 
handle 5-point vertex diagrams; (ii) proper account of summation over graviton excited modes; 
and (iii) modifications to handle massless spin-2 particles. As a result, processes with spin-2 
particles can now be simulated in a fully automatic way. 

Second, we have studied the simulation of graviton production in inclusive multi-jet samples. 
The fc^-MLM parton shower matching scheme has been employed for several benchmark models 
and points: for ADD we have chosen A = 5 TeV for the LHC, A = 1 TeV for the Tevatron, with 
(5 = 2, 4, 6; for the RS model, we have taken nii = 100 GeV and 1 TeV, and A = 3 TeV; and 
finally for MGM model, we defined M(/i^) = 1 and 2 TeV. We have considered pp/pp G + n- 
jets in all the scenarios. In ADD and MGM models, the graviton appears as missing energy while 
in RS model it may be observed through its decay products (here, example, we have 

chosen a pair of leptons). The fully inclusive samples have been normalized to the corresponding 
NLO inclusive total cross sections for producing G-\-X. Key distributions for the most relevant 
observables involving one or more jets have then been compared to those obtained by the NLO 
calculation G -|- jet -|- X. 

The overall agreeement in shape as well as in normalization between NLO observables and 
the corresponding ones obtained by MadGraph/MadEvent -|- Pythia is excellent. In par- 
ticular the fact that a unique normalization factor (the one normalizing the overall inclusive 
sample) suffices to describe not only the shape but also the normalization of jet distributions 
makes the use of matched samples in experimental analyses straightforward and accurate. Fi- 
nally, we note that observables which are sensitive to the number of jets, such as the Ht are 
better described by the matched sample than a fixed-order calculation. 

In conclusion we have implemented and validated the generation of fully inclusive samples 
with the correct leading-order jet multiplicity for spin-2 particle production at the Tevatron 
and the LHC. We look forward to their use in the current and forthcoming searches at hadron 
colliders. 
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A New HELAS subroutines for the 5-point vertex implementation 

Recall that ggggT vertices are obtained from the interaction Lagrangian among the graviton 
tensor and four gluons by Eq. (3.1). The subroutine GGGGTX computes the portion of the ampli- 
tude of the ggggT vertex from four Gluon polarization vectors and a Tensor boson wavefunction 
corresponding to the color structure J^^'eycde^ which should be called as 

CALL GGGGTX(VA, VB, VC, VD, TC, GC, GT , VERTEX) (A.l) 

The inputs VA(6) , VB(6) and VC(6) are complex six-dimensional arrays which contain the Vector 
boson wavefunctions, and their momenta as 

= (SReVA(5),KeVA(6),9?nVA(6),3=mVA(5)), 
= (3?eVB(5),5ReVB(6),3=mVB(6),9mVB(5)), 
= (3?eVC(5),KeVC(6),3=mVC(6),3=mVC(5)). 

The input TC(18) is a complex 18-dimensional array which contains the wavefunction of the 
Tensor boson: 

T^^+^'^+i =TC(4// + i/ + l), (A.2) 

and its four-momentum as 

p^" = (KeTC(l7),KeTC(18),9?nTC(17),9mTC(l7)), 
The output VERTEX is a complex number in units of GeV: 

VERTEX = -GT GC^ T"" G^^^^pXaS V^V^V^V^ (A.3) 

with 

Gpiy^paxs = VtiiyiVpcrVxs - VpSritTx) 

- ^pxVuSVpa + (pt ^ i^)] , (A.4) 

and we use the notation 

y^ = vA(/z + i), y^ = vB(^ + i), 

= VC(Ai + 1), Vi; = VD(Ai + 1). (A.5) 
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In order to insert automatically 5-point vertices in MadGraph, we introduce non-propagating 
colorless tensor boson tA as an auxiliary particle. With the portion of the 5-point vertex with 
color structure f^-^'^f'^'^'^ can be reduced automatically. After introducing particle Ia^ several new 
HELAS subroutines were added. They are summarized on table 1, and can be called on the 
following way: 

VVTCXX 

This subroutine computes the amplitude of the ggtA vertex, 

CALL VVTCXX(VA,VB,TC,GC,VMASS , VERTEX) 

VMASS represents the vector gluon mass which is zero (although it does not play any role here, 
we keep it as an input argument in accordance with the subroutine VVTXXX, for convenience of 
MadGraph). What we compute here is 

VERTEX = GC ri^pT^^^T^'^V^V^. (A.6) 

JVTCXX 

This subroutine computes an off-shell vector current J made from the interactions of a Vector 
gluon and an auxiliary Tensor boson by the ggtA vertex, and should be called as 

CALL JVTCXX(VC,TC,GC, VMASS, VWIDTH , JVTC) 

The input VC(6) is the wavefunction and momentum of the gluon. The output JVTC (6) gives 
the off-shell vector current multiplied by the gluon propagator, which is expressed as a complex 
six-dimensional array: 

JVTC(a + 1) = -,^„^y°r^-, (A.7) 

r 

and 

JVTC(5) = V(5) +TC(17), (A. 8) 

JVTC(6) = V(6) + TC(18). (A.9) 

Here the momenta q are 

= (5ReJVTC(5),3?eJVTC(6),9mJVTC(6),9mJVTC(5)). 

uvvcxx 

This subroutine computes an off-shell tensor current U for the auxiliary tensor tA, made from 
two flowing-out Vector gluons by the ggt^ vertex, and should be called as 

CALL UVVCXX(VA, VB, GC, VMASS, XM, XW , UVVC) 

The inputs XM and XW are two dummy arguments for which we keep them in accordance with 
the subroutine UVVXXX. The output UVVC (18) is a complex 18-dimensional array: 

T'^P = GCV^V^ (A.IO) 
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for the first 16 components of UVVC, and 

UVVC(17) = VA(5) + VB(5), (A. 11) 

UVVC(18) = VA(6) + VB(6). (A. 12) 

TTTXXX 

This subroutine computes the ampHtude of the tA^A^ vertex, called as 

CALL TTTXXX(TC,T1C,T2C,GT, VERTEX) 

The inputs TIC (18) and T2C(18) are complex 18-dimensional arrays which contain the wave- 
function and momenta for the auxiliary tensors. What we compute here is 

VERTEX = GT T^^T^'^TfGf.^^p^xs, (A. 13) 

with Ti^2 defined from TIC (18) and T2C(18) as in Eq. (A. 2). 

UTTAXX 

This subroutine computes an off-shell non-propagating tensor current U, made from the flowing- 
out graviton Tensor and the auxiliary Tensor by the tA^AT vertex, and should be called as 

CALL UTTAXX(TC, TIC, GT, UTTA) 
The output UTTA (18) is a complex 18-dimensional array: 

Tf = GT r?-"r?'^''r^^rf^G^,,p„A5 (A. 14) 

for the first 16 components of UTTA, and 

UTTA(17) = T(17) -FT1(17), (A. 15) 

UTTA(18) = T(18) -FT1(18). (A.16) 

UTTBXX 

This subroutine computes an off-shell graviton tensor current U, made from the two flowing-out 
auxiliary Tensors by the t^T vertex, and should be called as 

CALL UTTBXX(T1C, T2C, GT, TMASS, TWIDTH, UTTB) 

The inputs TMASS and TWIDTH are the graviton mass and width, uit and Tt- The output 
UTTA (18) is a complex 18-dimensional array: 

T^^ = GT ^ "f^"'':^ ^ T,P^TfG,,,,^xs (A.17) 

for the first 16 components of UTTB, and 

UTTB(17) = T1(17) -FT2(17), (A.18) 
UTTB(18) = T1(17) -FT2(18). (A.19) 

Here q is the momentum of the off-shell tensor boson given in Eqs. (A.18) and (A.19) as 

q^" = (3f?eUTTB(l7),3fieUTTB(l8),9mUTTB(l8),3=mUTTB(l7)). 

And B^a,vp is the polarization summation tensor [16]. 



- 25 - 



B New subroutines for summation over continuous mass spectrum 

In order to have the mass integration reahzed in MadGraph within the framework of the ADD 
model, new auxihary particles were introduced in Sec. 3. In consequence, three new sub-routines 
summarized in table 2 were added into the HELAS structure. The details for the insertion of 
these sub-routines are the following: 

PXXXXX 

This subroutine stores the helicity and momentum of the auxiliary spin-2 Pseudo-particle for 
further usage, and should be called as 

CALL PXXXXX(P, XM, NHEL, NST , PC) 

The input P(0:3) is a real four-dimensional array which contains the four-momentum of the 
pseudo particle, NHEL (= ±2, ±1,0) specifies its helicity A, NST specifies whether the boson is in 
the final state (NST = 1) or in the initial state (NST = -1). XM is its mass but does not play any 
role here. 

The output PC (18) is a complex 18-dimensional array, among which only the following 
matter, namely 

PC( 1) = NHEL, (B.l) 

and 

(PC(17), PC(18)) = NST (P(0) + iP(3), P(l) + iP(2)). (B.2) 

TPSXXX 

This subroutine computes the amplitude of the X1X2T vertex, and should be called as 

CALL TPSXXX(TCO,PC,SC,GP,TMASSO,TWIDTH, VERTEX) 

SC(3) is a complex three-dimensional array which contains the momentum for xi. GP(1:2) here 
are not coupling constant, but contain the ADD model inputs: 

GP( 1) = A + i(5, 

GP( 2) = Mio^ + iMup, (B.3) 

Mioyj and M^p are the lower and upper limits for graviton mass integration. TMASSO (rriTo for 
below) is only temporarily used and got from the MG input file param_card.dat, while the real 
graviton mass ttit is defined as the invariant mass of xl and x2, i.e., 

mT^\qi + q2\, (B.4) 
with qi and (72 are the momenta of xi and X2. 
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Also note the input TCO is indeed calculated from the other side of Feynman diagram, and 
contains the polarization summation tensor B^yp^, which can be used to project out the real 
graviton wave function, due to the relations of the 5 helicity states of the tensor boson [16]: 

B^u,af}ip)= X] '^t^y{P^^)'^ap{p,^T, (B.5) 

A=±2,±l,0 

and 

6^^(p,A)e^,(p,A')* = 5AA', (B.6) 

with p = qi + q2, and the helicity A is got from the pseudo-particle's wavefunction (Eq. (B.l)). 
Now we can get the true graviton tensor wave function by calling the subroutine 

CALL TXXXXX(P, mT, A, +1, TC) (B.7) 

Finally, the output VERTEX is 

VERTEX =^^p{mT){p'^ - + imToTT)TavTlt'' 
rriT 

xe{mT- Mion,)0iMup-mT), (B.8) 

where 9 represents the Heaviside step function. In Eq. (B.8), we include the mass density 
factor(Eq. (3.2)), the inverse graviton propagator, and the compensation factor for the decay 
phase space. 

UPSXXX 

This subroutine computes an off-shell graviton tensor current U, made from the flowing-out 
auxiliary Scalar and Pesudo-particles, by the X1X2T vertex, and should be called as 



CALL UPSXXX(PC, SC, GP, TMASSO, TWIDTH, UPS) 

The output UPS is a complex 18-dimensional array, which is indeed got simply by calling inside 
UPSXXX 

CALL TXXXXX(P, mr, A, +1, UPS). (B.9) 
and including the step functions in Eq. (B.8). 



C New subroutines for the massless spin-2 particles 

In order to have massless spin-2 particles also implemented in MadGraph/MadEvent , the 
tensor wave function subroutine was modified, and four subroutines for the off-shell tensor were 
added to the former HELAS implementation. The details are presented as follows. 
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TXXXXX 

This subroutine computes the graviton tensor wave function, namely e^'^{p,X) and e^^{p,\)* in 
terms of the graviton four-momentum p and its hehcity A. The sum of hehcity states for the 
graviton is introduced on the routines to be given by Eq. (B.5). While the massive graviton 
has five helicity states of polarization, the massless spin-2 particles will only have two physical 
states. The sum of polarization for the massless graviton is therefore modified on the routines 
to set A = ±1, to zero allowing only A = ±2 to contribute. 

UIOXXX, UVVXXX, UIOVXX, uvvvxx 

These subroutines compute off-shell tensor currents U by the FFT, VVT, FFVT, and VVVT vertices 
respectively. The main modification here was the inclusion of the zero mass graviton propagator 
to be the one showed in Eq. (3.4) for the massless case. 

D Higgs effective theory new implementation into MadGraph 

In Ref. [15], Higgs effective theory has been implemented into MadGraph by reorganizing the 
Lagrangian with introducing an extra non-propagating auxiliary tensor, to avoid the 5-point 
vertex problem. Although it seems more elegant in physics, it is hard to find a similar way for 
graviton case (due to more complicated Lagrangian for graviton interactions). Moreover, the 
previous way affects the 4-gluon vertices, reducing them into subparts corresponding to different 
color structure, and thus increases the numbers of Feynman diagrams for process like gg — )• gg. 
This may slow the running of MadEvent, which is based on single diagram enhanced method 
and thus sensitive to the singular diagram numbers. 

As mentioned in Sec. 3, we have implemented HEFT into MadGraph in the same way as 
for the graviton (see the model directory heftb in MadGraph). We add three more HELAS 
subroutines as following: 

TTSCXX 

This subroutine computes the amplitude of the Ia^a^ vertex, 

CALL TTSCXX(T1C, T2C, S, GH, VERTEX) 

The input GH = as/{'iiTv) is the Higgs effective coupling constant. S(3) is a complex three- 
dimensional array which contain the wavefunctions of the Scalar bosons, S(l), and their four- 
momenta as 

p^ = (KeS(2),KeS(3),$>mS(3),3mS(2)). 

The output VERTEX is: 

VERTEX = -GH5(1) T^^Tf {i^^pr^^, - 7?^„r?,^) (D.l) 
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HTTCXX 

This subroutine computes an off-shell scalar current H made from the interactions of two auxiliary 
Tensor bosons by the tA^AH vertex, and should be called as 

CALL HTTCXX(T1C, T2C, GH, SMASS, SWIDTH, HTTC) 

The inputs SMASS and SWIDTH are the scalar boson mass ms and width Ts- The output HTTC (3) 
gives the off-shell scalar current multiplied by the scalar boson propagator, which is expressed 
as a complex three-dimensional array: 

HTTC(l) = — —T^^Tfir^^,^^^ - r],^) (D.2) 

— rrig + i7nsTs 

and 

HTTC(2) = T1C(17) + T2C(17), (D.3) 
HTTC(3) = T1C(18) +T2C(18). (D.4) 

Here the momenta q is 

= (KeHTTC(2),3fteHTTC(3),9mHTTC(3),3=mHTTC(2)). 

UTSCXX 

This subroutine computes an off-shell auxiliary tensor current U made from one flowing-out 
Scalar bosons and one auxiliary Tensor boson by the t^tXii vertex, and should be called as 

CALL UTSCXX(T1C, S, GH, XM, XW , UTSC) 

The inputs XM and XW are two dummy arguments. The output UTSC (18) gives the off-shell tensor 
current, which is expressed as a complex 18-dimensional array: 

T''^ = -GH5(1) rr(r/^/3r?a. - V,.aV.p) (D-5) 

for the first 16 component of UTSC, and 

UTSC(17) = T1C(17) S(2), (D.6) 
UTSC(18) = T1C(18) S(3). (D.7) 



- 29 - 



References 

[1] I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos, and G. R. Dvali, New dimensions at a millimeter 
to a Fermi and superstrings at a TeV, Phys. Lett. B436 (1998) 257-263, [hep-ph/9804398]. 

[2] N. Arkani-Hamed, S. Dimopoulos, and G. R. Dvali, Tlie tiierarchy problem and new dimensions at 
a millimeter, Phys. Lett. B429 (1998) 263-272, [hep-ph/9803315]. 

[3] N. Arkani-Hamed, S. Dimopoulos, and G. R. Dvali, Phenomenology, astrophysics, and cosmology 
of theories with submillimeter dimensions and TeV scale quantum gravity, Phys. Rev. D59 (1999) 
086004, [hep-ph/9807344]. 

[4] G. Dvali, G. Gabadadzc, M. Kolanovic, and F. Nitti, Scales of gravity, Phys. Rev. D65 (2002) 
024031, [hep-th/0106058]. 

[5] G. Dvali, Black Holes and Large N Species Solution to the Hierarchy Problem, Fortsch. Phys. 58 
(2010) 528-536, [arXiv: 0706 . 2050]. 

[6] G. Dvali and M. Rcdi, Black Hole Bound on the Number of Species and Quantum Gravity at LHC, 
Phys. Rev. D77 (2008) 045027, [arXiv:0710.4344]. 

[7] X. Calmet and S. D. H. Hsu, TeV gravity in four dimensions?, Phys. Lett. B663 (2008) 95-98, 
[arXiv: 0711. 2306]. 

[8] X. Calmet, S. D. H. Hsu, and D. Reeb, Quantum gravity at a TeV and the renormalization of 
Newton's constant, Phys. Rev. D77 (2008) 125015, [arXiv : 0803 . 1836]. 

[9] X. Calmet, S. D. H. Hsu, and D. Recb, Grand unification and enhanced quantum gravitational 
effects, Phys. Rev. Lett. 101 (2008) 171802, [arXiv : 0805 . 0145]. 

[10] L. Randall and R. Sundrum, A large mass hierarchy from a small extra dimension, Phys. Rev. 
Lett. 83 (1999) 3370-3373, [hep-ph/9905221]. 

[11] L. Randall and R. Sundrum, An alternative to compactification, Phys. Rev. Lett. 83 (1999) 
4690-4693, [liep-th/9906064]. 

[12] J. M. Buttcrworth et. at. The Tools and Monte Garlo working group Summary Report, 
arXiv: 1003. 1643. 

[13] T. Stelzer and W. F. Long, Automatic generation of tree level helicity amplitudes, Gomput. Phys. 
Commun. 81 (1994) 357-371, [hep-ph/9401258]. 

[14] F. Maltoni and T. Stelzer, MadEvent: Automatic event generation with MadGraph, JHEP 02 
(2003) 027, [hep-ph/0208156]. 

[15] J. Alwall et. al, MadGraph/MadEvent v4: The New Web Generation, .JHEP 09 (2007) 028, 
[arXiv: 0706. 2334]. 

[16] K. Hagiwara, J. Kanzaki, Q. Li, and K. Mawatari, HELAS and MadGraph/MadEvent with spin-2 
particles, Eur. Phys. J. C56 (2008) 435-447, [arXiv : 0805 . 2554]. 

[17] S. Karg, M. Kramer, Q. Li, and D. Zeppenfeld, NLO QCD corrections to graviton production at 
hadron colliders, Phys. Rev. D81 (2010) 094036, [arXiv : 0911 . 5095]. 

[18] P. Mathews, V. Ravindran, and K. Sridhar, NLO-QCD Corrections to Dilepton Production in the 
Randall- Sundrum Model, JHEP 0510 (2005) 031, [hep-ph/0506158]. 

[19] Q. Li, C. S. Li, and L. L. Yang, Soft Gluon Resummation Effects in Single Graviton Production at 
the CERN Large Hadron Gollider in the Randall- Sundrum Model, Phys. Rev. D74 (2006) 056002, 
[hep-ph/0606045]. 



- 30 - 



S. Catani, F. Krauss, R. Kuhn, and B. Webber, QCD Matrix Elements + Parton Showers, JHEP 
0111 (2001) 063, [hep-ph/0 109231]. 

F. Krauss, Matrix elements and parton showers in hadronic interactions, JHEP 0208 (2002) 015, 
[hep-ph/0205283]. 

S. Hoeche et. ai, Matching parton showers and matrix elements, hep-ph/0602031. 

J. Alwall et. ai, Comparative study of various algorithms for the merging of parton showers and 
matrix elements in hadronic collisions, Eur. Phys. J. C53 (2008) 473-500, [arXiv: 0706. 2569]. 

J. Alwall, S. de Visscher, and F. Maltoni, QCD radiation in the production of heavy colored 
particles at the LHC, JHEP 02 (2009) 017, [arXiv : 0810 . 5350]. 

J. Campbell, Next-to-leading order corrections to W-l-2 jet and Z+2 jet production at hadron 
colliders, Phys. Rev. D65 (2002) 113007. 

T. Han, J. D. Lykken, and R.-J. Zhang, On kaluza-klein states from large extra dimensions, Phys. 
Rev. D59 (1999) 105006, [liep-ph/98 11350]. 

G. F. Giudice, R. Rattazzi, and J. D. Wells, Quantum gravity and extra dimensions at high-energy 
colliders, Nucl. Phys. B544 (1999) 3-38, [hep-pli/9811291]. 

E. A. Mirabelli, M. Perelstein, and M. E. Peskin, Collider signatures of new large space 
dimensions, Phys. Rev. Lett. 82 (1999) 2236-2239, [hep-ph/9811337]. 

D. Kapner, T. Cook, E. Adelberger, J. Gundlach, B. Heckel, C. Hoyle, and H. Swanson, Tests of 
the gravitational inverse- square law below the dark-energy length scale, Phys. Rev. Lett. 98 (2007) 
021101, [hep-ph/0611184]. 

N. Kaloper, J. March-Russell, G. D. Starkman, and M. Trodden, Compact hyperbolic extra 
dimensions: Branes, kaluza-klein modes and cosmology, Phys. Rev. Lett. 85 (2000) 928-931, 
[hep-ph/0002001]. 

K. R. Dienes and A. Mali, Shadows of the Planck scale: The changing face of compactification 
geometry, Phys. Rev. Lett. 88 (2002) 111602, [hep-th/0111264]. 

G. F. Giudice, T. Plehn, and A. Strumia, Craviton collider effects in one and more large extra 
dimensions, Nucl. Phys. B 706 (2005) 455-483, [hep-ph/0408320]. 

LEP Exotica Working Group, LEP Exotica WC 2004-03. 

CDF Collaboration, T. Aaltonen et. ai. Search for large extra dimensions in final states 
containing one photon or jet and large missing transverse energy produced in pp collisions at ^/s = 
1.96-TeV, Phys. Rev. Lett. 101 (2008) 181602, [arXiv : 0807 . 3132]. 

DO Collaboration, V. M. Abazov et. ai. Search for large extra dimensions via single photon plus 
missing energy final states at y/s ~ 1.96-TeV, Phys. Rev. Lett. 101 (2008) 011601, 
[arXiv: 0803. 2137]. 

CDF Collaboration, Search for Large Extra Dimensions in the Production of Jets and Missing 
Transverse Energy in pp Collisions at ^/s — L96 TeV, Phys. Rev. Lett. 97 (2006) 171802, 
[hep-ex/0605101]. 

H. Davoudiasl, Echoes from a Warped Dimension, Nucl. Phys. Proc. Suppl. 200-202 (2010) 
149-158, [arXiv: 0909. 1587]. 

CDF Collaboration, T. Aaltonen et. ai, A search for high-mass resonances decaying to dimuons 
at CDF, Phys. Rev. Lett. 102 (2009) 091805, [arXiv : 0811 . 0053]. 



- 31 - 



[39] DO Collaboration, V. M. Abazov et. al, Search for Randall- Sundrum gravitons with 1 fb ^ of data 
frompp collisions at yfs = 1.96- TeV, Phys. Rev. Lett. 100 (2008) 091802, [arXiv : 0710 . 3338]. 

[40] F. Larsen and F. Wilczek, Renormalization of black hole entropy and of the gravitational coupling 
constant, Nucl. Phys. B458 (1996) 249-266, [hep-tli/9506066]. 

[41] X. Calmet and P. de Aquino, Quantum Gravity at the LHC, Eur. Phys. J. C68 (2010) 305-311, 
[arXiv: 0906. 0363]. 

[42] X. Calmct and M. Feliciangeli, Bound on four- dimensional Planck mass, Phys. Rev. D78 (2008) 
067702, [arXiv: 0806. 4304]. 

[43] X. Calmct, P. dc Aquino, and T. G. Rizzo, Massless versus Kaluza-Klein Gravitons at the LHG, 
Phys. Lett. B682 (2010) 446-449, [arXiv : 0910 . 1535]. 

[44] H. Murayama, I. Watanabc, and K. Hagiwara, "HELAS: HELicity amplitude subroutines for 
Feynman diagram evaluations." KEK-91-11, 1991. 

[45] K. Hagiwara, H. Murayama, and I. Watanabe, Search for the Yukawa interaction in the process e-h 
e- -> t anti-t Z at TeV linear colliders, Nucl. Phys. B367 (1991) 257-286. 

[46] K. Hagiwara, P. Konar, Q. Li, K. Mawatari, and D. Zeppenfeld, Graviton production with 2 jets at 
the LHC in large extra dimensions, JHEP 04 (2008) 019, [arXiv : 0801 . 1794]. 

[47] H. van Dam and M. J. G. Veltman, Massive and massless Yang-Mills and gravitational fields, 
Nucl. Phys. B22 (1970) 397-411. 

[48] T. Sjostrand, S. Mrenna, and P. Skands, Pythia 6.4 physics and manual, JHEP 0605 (2006) 026, 
[hep-ph/0603175]. 

[49] G.Corcella, I.G.Knowles, G.Marchesini, S.Moretti, K.Odagiri, P.Richardson, M.H.Seymour, and 

B.R. Webber, Herwig 6.5: an event generator for hadron emission reactions with interfering gluons 
(including supersymmetric processes), JHEP 0101 (2001) 010, [hep-ph/00 11363]. 

[50] F. Krauss, A. Schaelicke, S. Schumann, and G. SofF, Simulating w/z-hjets production at the 
tevatron, Phys. Rev. D70 (2004) 114009, [hep-ph/0409106]. 

[51] S. Mrenna and P. Richardson, Matching matrix elements and parton .showers with herwig and 
pythia, JHEP 0405 (2004) 040, [hep-ph/0312274]. 

[52] L. Vacavant and I. HinchlifFe, Signals of models with large extra dimensions in ATLAS, J. Phys. 
G27 (2001) 1839-1850. 

[53] J. Pumplin et. al.. New generation of parton distributions with uncertainties from global QCD 
analysis, JHEP 07 (2002) 012, [hep-ph/0201195]. 

[54] M. Rubin, G. P. Salam, and S. Sapeta, Giant QGD K-factors beyond NLO, JHEP 09 (2010) 084, 
[arXiv: 1006. 2144]. 



- 32 - 



